golv-int/ 9706005 



A Lax Description for Polytropic Gas Dynamics 



J. C. Brunelli* 

Universidade Federal de Santa Catarina 
Departamento de Fisica - CFM 
Campus Universitario - Trindade 
CP. 476, CEP 88040-900 
Florianopolis, SC - BRAZIL 

and 

Ashok Das 

Department of Physics and Astronomy 
University of Rochester 
Rochester, NY 14627 - USA 



Abstract 

We give a Lax description for the system of polytropic gas equations. The special 
structure of the Lax function naturally leads to the two infinite sets of conserved charges 
associated with this system. We obtain closed form expressions for the conserved charges 
as well as the generating functions for them. We show how the study of these generating 
functions can naturally lead to the recursion relation between the conserved quantities as 
well as the higher order Hamiltonian structures. 
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Systems of hydrodynamic type have been extensively studied over the last several 
years [1-5]. These are first-order quasilinear systems of partial differential equations. A 
sub class of these systems is the two-component hyperbolic system of conservation laws [5] 
which has the general Hamiltonian form (in 1 + 1 dimensions) 
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is an Hamiltonian functional of zeroth order. An infinite number of zeroth order conserved 
charges are associated with this system of equations. In fact, if 



Q[it,u] = J dxQ(u,v) 



(4) 



represents a conserved charge of the system, it is easy to show that the density must satisfy 

2 H 8 2 Q _ 8 2 H 8 2 Q 
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In particular, we note that H satisfies this equation. 

Of special interest are systems for which (5) admits a separation of variables such that 

d 2 Q X(u) d 2 Q 



du 2 jj,(v) dv 2 



A separable Hamiltonian of the form 



1 



H(u,v) = - ( -u 2 v + F(v) 



(6) 



(7) 



describes the physically important system of gas dynamics [6] corresponding to A (it) = 1 
and fj,(v) = F"{y)/v. There is yet another important system - the nonlinear elastic medium 
model [6] - which is also described by a separable Hamiltonian of the form 

T 



H(u,v) = - -u 2 + F(v) 



(8) 
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Olver and Nutku [5], in a beautiful paper, gave a systematic description of systems 
like (1). In particular, they have discussed the conserved charges and the Hamiltonian 
structures of the system of polytropic gas, which corresponds to the choice 

F(v) = - { ^- Ty 7 ^0,1 (9) 

in (7), leading to the equations 

U t + UU X + V^^Vx =0 

(10) 

vt + (uv) x =0 

It was shown in [5] that this system of equations has two infinite sets of conserved charges. 

However, a Lax description for such systems is missing so far. It is known [1,3] that a 
Lax description for dispersionless systems such as (1) are best given in terms of polynomial 
functions in phase space (We will continue to call these as Lax operators.). However, 
whatever Lax formalism describes such a system, it must have the peculiar property of 
yielding two infinite sets of conserved charges. In Ref. 7 a nonstandard Lax description for 
this system of equations with 7 = 2 was proposed, which is also equivalent to the shallow 
water wave equation. This system of equations is also known as the irrotational Benney's 
equation [4] and have been well studied in the literature. However, 7 = 2 turns out to be 
a particularly simple choice in that the two sets of conserved charges become degenerate 
in this case and hence one does not expect any new structure in the Lax description. In 
going beyond 7 = 2 to 7 = 3, one immediately notices that the system of equations (10) 
can be written as the 2x2 matrix Riemann equation 

U t = UU X (11) 

where 

V=( U V ) (12) 

It is quite straightforward to show that the dispersionless Lax operator L = p 2 l + U, 
describes this system of equations with the two infinite sets of conserved charges obtained 
from TrL - ^ = j dxRes (trL - ^ j taking the matrix trace "tr" both with respect to 



diagonal and off-diagonal elements of the matrices ("Res" is the coefficient of the p~ x 
term). 

For any integer 7, therefore, either a matrix or a scalar generalization of a Lax repre- 
sentation for (10) can be pursued. We have found that the scalar Lax operator 

L = p 7-l +w+ _^L_ p -(7-l) (13) 

with the dispersionless nonstandard representation 

where {A, B} = — §f leads to the system of polytropic gas equations (10) (In 

fact, the hierarchy of equations is obtained from ). It is interesting to note 

that the Lax operator (13) as well as the dynamical equations (14) naturally reduce to the 
shallow water system for 7 = 2 given in Ref. [7]. However, we would like to emphasize 
here that the Lax function in (13) is quite rigid in that any deformation of this leads to 
inconsistent equations. It is also worth noting here that since the system of equations 
(10) is known to have two distinct recursion operators [5], we suspect that there may be 
yet another equivalent Lax description of the system. However, we have not succeeded in 
finding it. 

Conserved charges associated with the system can be obtained from (13) through 
(These are the only fractional powers with nontrivial residues.) 

U n = Tr L n+ ^r , n = 0, 1, 2, 3, . . . (15) 

and the first few densities are easily obtained to be 
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H 2 = (37-4)(27- 3)(7-2) / I 1 , 

(7-1)3 ^3, 7 _! 7 _ 2 
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In general, we can write 

Q + l)(7-l)-l 

H n = (n + l)\C n+1 ^- 1) H. 
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where H^s correspond to one of the infinite sets of conserved densities obtained in [5]. The 
conserved densities (16) were obtained by expanding L^ 1 around p = oo. However, the 
present Lax function has a singularity at p = and, consequently an alternate expansion 
of the fractional powers is possible around p = as well [8] and it gives us a second set of 
nontrivial conserved charges through 

H n = TrL n+ ^, n = 0,1,2,3,... (17) 

where the first few densities are 



Hi =( 7 -l)-^T ^^ m; 
5>=(7-l)-*fe)£(^. + 
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In general, we can write 



TT p (T-l) TT 

n — 2 11 n 

(7 - 1) -y- 1 



where H^s correspond to the second set of conserved densities obtained in [5]. Incidentally, 
this alternate expansion around p = also gives us a second consistent dispersionless Lax 
equation 

% = Ul&)..,l) (19) 



dt [\ ) <o 
which leads to the set of equations 

u t = - ' 



.(7"2)Vx (20) 
With proper rescaling, this can be rewritten also as 

T-3 

Ut = - V 7 V x 



We recognize this as the elastic medium equations with the polytropic gas potential F(v) = 

7 — 1 

( 7 J"i)( 7 _ 2 ) [5] • And our derivation clarifies the fact that both these systems share the same 
conserved charges because they are obtained from the same Lax function. 

We see that the Lax operator (13) nicely reproduces both sets of densities H n and H n 
described by Olver and Nutku [5] . We can even obtain a closed form expression for these 
densities from their definition and they have the form 
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Since these are obtained from a Lax function, it follows that they are conserved under the 
flow. However, it is also straightforward to check explicitly from the closed form of these 
charges in (21) that 



d 2 H n 3 d 2 H n 
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dv 2 du 2 ' dv 2 du 2 

leading (from (5)) to the fact that they are indeed conserved. 

These expressions allow us to write the following generating functions [1,3,4] for the 
conserved densities (In fact, they also follow from an analysis of the associated linear 
equations in the dispersionless limit.) 
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We note that these generating functions reduce to those of Manin [1] for 7 = 2. Further- 
more, using (10), it is easy to check that 
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This shows that x an d X are ? indeed, conserved for any value of the parameter A and, 
consequently, generate conserved densities. In fact, expanding these for large A, it is easy 
to identify 

00 C_1 ) n \-( n + 1 ) ("+1)(t -1)-1 

n-0 , 2A s 

^(n( 7 -l) + l) ( 7 -l) — 
It is interesting to note that the roots of the radical in (22) give the Riemann invariants 
of the problem, namely, 

7-1 

x 2^ 

—A± = U ± 



7-1 

and that both the polytropic gas and the nonlinear elastic medium with a polytropic gas 
potential share the same Riemann invariants. 

The generating functions satisfy various identities. It is easy to check from the defi- 
nitions in (22) that Q = 
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It follows trivially from these that 
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The consistency of these leads to 

^ =W T-3^X 5 ^X =w 7-3^ (27) 

<9i> 2 <9-u 2 ' <9i> 2 du 2 

which is indeed in the form (6). We can now prove that the conserved charges are in 
involution with respect to the Hamiltonian structure in (1) in a simple manner [4]. By 
definition 



{x(a),x(a )} =j d X + -^r^^r 

= J dx (Fu x +Gv x ) 



(28) 



where 

F ^dxW d 2 X {\>) + dxW d 2 X {\') 
du dudv dv du 2 
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Using (27) we can show that 
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which implies that 



Consequently, we obtain 



F= M£ G= dK 
du ' dv 



{X (A), X (A<)} = fd* (f «. + f «.) = /dxf = (31) 
Similarly, one can show 

{X(A),X(A / )} = {X(A),X(A')} = (32) 

Since the generating functions are in involution for arbitrary parameters A and A', it follows 
that all the conserved charges, H n and Tt n , must also be in involution with one another. 

The generating functions, x an d X, contain all the information about the system. One 
can easily study the symmetries, the recursion relation between the conserved quantities as 
well as the higher order Hamiltonian structures from these. Here we give the simple, but 
nontrivial example of how the recursion relations between the conserved charges as well 
as the second Hamiltonian structure for the system can be derived from the generating 
function. Let 

1- = ^EA-"!- 03) 

Then, we can write the entire hierarchy of polytropic gas equations in the compact form 
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From Eq. (25), we can write 



(7-l)(« + A)^ + 2t;^-x = 

9V 8 ~ (35) 
(7-1) (2 V -(. + A))g-2^- 2 g = 
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Taking the derivative of these with respect to x, we obtain 

-(7 - l)Xd x ^ =(7 - 2)u x ^ +v x ^ + ( 1 - l)ud x ^- + 2v d x ^- 
w ; du w ; dw w ' du dv ^ 

-(7 - 1)A0*|* =(7 - 2 K~ 3 ^ + 2vl ~ 2d -^ + u -% + (7 - ^"f* 

Restricting to specific powers of A, these, then define a recursion relation between the 
conserved densities (with suitable normalization), namely, 
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We note here that the recursion relation in (36) can equivalently be written as 
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where 



( dv^~ 2 + vi- 2 d du+( 1 - 2)ud \ 
1 \{ 1 -2)du + ud dv+vd J [ } 



Using this in the dynamical equation, we note that we can write 

A" 1 






(40) 



which shows that the hierarchy of equations is bi-Hamiltonian with T>\ representing the 
second Hamiltonian structure obtained in [5] . 

Finally, we note here that the system of polytropic gas equations has a natural zero 
curvature representation in this framework [3]. Let 

then, it is easy to check that 

d t Ai-d x A o -{A o ,Ai} = (42) 



gives the equations for the polytropic gas in (10). 
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